Double Inflation and the Low CMB Quadrupole 
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Recent released WMAP data show a low value of quadrupole in the CMB temperature fluctua- 
tions, which confirms the early observations by COBE. In this paper we consider a model of two 
inflatons with different masses, V((fri,(f>2) = |mf</>f + \m\(fy , mi > 7712 and study its effects on 
CMB of suppressing the primordial power spectrum P(k) at small k. Inflation is driven in this model 
firstly by the heavier inflaton cj>\, then the lighter field <j>2- But there is no interruption in between. 
We numerically calculate the scalar and tensor power spectra with mode by mode integrations, then 
fit the model to WMAP temperature correlations TT and the TE temperature-polarization spectra. 
Our results show that with mi ~ I0 14 GeV and ni2 ~ 10 13 GeV, this model solves the problems of 
flatness etc. and the CMB quadrupole predicted can be much lower than the standard power-law 
ACDM model. 
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Recently the Wilkinson Microwave Anisotropy Probe 
(WMAP) data [J H have been released and 
it is shown that the data is consistent with the predic- 
tions of the standard ACDM model with an almost scale- 
invariant, adiabatic and Gaussian primordial (scalar) 
fluctuations. However, there remain intriguing discrep- 
ancies between the model and the observations, which 
show the overprediction of the model on the ampli- 
tudes of fluctuations at both the largest and the small- 
est scales. In Ref.0 Spergel et al. include other data 
of the Cosmic Microwave Background(CMB)0, 0] and 
Large Scale Structure(LSS)[l They find that for 
power law ACDM model the best fit for the amplitude 
of fluctuations gradually drops as the probe of scale k 
increases and the data supports for a nonzero running of 
the scalar spectrum index from blue to red at 2a with 
dn s /dlnk = -0.042±g;8|J 0. In Ref.[I3|, the authors 
have questioned about the validity of the use of Lyman- 
a forest data despite this a slightly running of the 
power spectrum index is still favored in the analysis of 
Refs. jl [Hill EE!. Especially, with a detailed recon- 
struction of the power spectrum Mukherjee and Wang[l2j 
have shown a preferred feature at k ~ 0.01 Mpc -1 , con- 
sistent with a running of the index. And in Ref. 11] Bridle 
et al. have stressed the importance of the data for the 
first three multipoles I = 2, 3,4 on the requirement for 
the running index. 

Theoretically there have been studies in the literature 
since the release of the WMAP data on models of in- 
flation whic h p rovide a running index required by the 
WMAP mmmmill]- However, as shown in Ref.0, 
the probability of finding a lower value of quadruple in 
the presence of a constant running of the spectral in- 
dex is no more than 0.9 percent for a spatial flat ACDM 
cosmology. Indeed, the lack of CMB power on the large 
angular scales seen already in COBE J2jj and rein- 



forced by WMAP is more challenging. This discrepancy 
may be due to cosmic variance. On the other hand, it 
probably gives a hint for new physics. 

Recently several possibilities to alleviate the low- 
multipoles pro blem have be en p roposed in the literature 
[immiimmEimi which include consider- 
ing the effect of a finite universe Q, the late time inte- 
grated Sachs- Wolf (ISW) effect by quintessence pi |24|. 
a non-flat universe [22], a suppression of the primor- 
dial fluctuations [Til l24l |2f| or different release of the 
d atapll |2f1 | . In the framework of inflation, Contaldi et al 
in|24j have discussed two approaches to suppressing the 
large scale power. Beside the one of changing the infla- 
ton potential, they have proposed another one of chang- 
ing the initial conditions at the onset of inflation relative 
to the standard chaotic inflation model j2?|]. For the lat- 
ter case, the inflaton has to be assumed in the kinetic 
dominated regime initially. 

In this paper we consider a double inflation model and 
study the possibility of suppressing the lower multiples 
in the CMB. For a quantitative investigation we study a 
modelHl: 
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The double inflation in the literature has been studied 
widely. And phenomenologically the model in (1) can 
be realized naturally in particle physics. For example, 
in the sneutrino inflation models [29l l30j|. there are three 
sneutrinos which belong to three different families. Tak- 
ing two of them degenerated, it is effectively a model of 
double inflation. 

We assume in (1) that <j>\ is heavier than 02, i-e. mi > 
TO2- The inflation is firstly driven by <pi, then by <f>2, and 
there is no interruption in between. The transition takes 
place at m ( f >1 ~ H, where H is the Hubble rate. When 
the transition happens <f>\ starts to oscillate around the 
minimum of its potential. We denote the wavenumbcr 
of comoving mode which crosses the horizon around this 
moment as kf. Choosing the model parameters so that 
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kf corresponds to the scale around our current horizon, 
we will show in this paper that model (1) provides a scalar 
power spectrum much suppressed around kf. With a set 
of the model parameters we will give a specific example 
of the initial power spectra and fit the spectra to WMAP 
data. Our results show that the spectrum with a feature 
is favored and lower CMB multipolcs can be achieved by 
the spectrum with a feature provided by this model. 

For the discussions on double inflation, we use the 
notations of Ref . [3ll| . In a spatially flat Friedmann- 
Robertson- Walker (FRW) universe the evolution of the 
background fields for the potential given in (1) is de- 
scribed by the Klein-Gordon equation: 



0, + 3H4>i + y 0J = o , 

and the Friedmann equation: 
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where / = 1,2, a is the scale factor, the dot stands for 
time derivative and V x = dV/dx. Scalar linear pertur- 
bations to the FRW metric can be expressed generally as 
(we use the metric convention +, — , — , — ): 

ds 2 = (l+2A)dt 2 ~2aB l dx l dt-a 2 [(l-2ijj)5 tj +2E tJ ]dx l dx 

Thus the equation for the evolution of the perturbation 
5<pi with comoving wavenumber k is given by 



3H 5(j), 
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Defining the adiabatic field a and its perturbation as [3jJ : 

(6) 



a = (cos 0)0i + (sin 0)02 , 
da = (cos 0)<50i + (sin 0)502 



with 1 



sin0 



n -r<n v V1 

The background equations (3) and (2) become 
8ttG ,1 
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H 



3<2 &2 + V ^ 
a + 3Ha + V a = , 



(8) 



where V a = (cos 0)1^ + (sin0)V^, . The comoving cur- 
vature perturbation is given by 



TZ = tp-\ 5a . 

a 



(9) 



1 Our definitions of cos 6 and sin 9 have the opposite signs with 
those in Ref. l3ll. but these differences do not affect the results. 



We assume that there is no entropy perturbation, and 
this is consistent with the results of WMAP .4] . Under 
this assumption, there is an adiabatic condition between 
501 and S(f>2- 



(10) 



So, the equation governing the evolution of adiabatic per- 
turbation is the same as that in the single field inflation 
model 0,113: 



(k 2 )u = 

z 



(11) 



where u = —zR, and z = aa/H, the prime denotes the 
derivative with respect to conformal time r\ (r\ = J — ). 
The power spectrum of adiabatic perturbation is defined 

as 



Pn(k) 



2tt z z 



(12) 



which approaches a constant at late time as k/aH — » 0. 
Similarly for tensor perturbations the power spectrum is 



Pg(k) 



k 3 v 



(13) 



where v — a '5 , is the linear tensor perturbation 3 
and the equation of motion for v is : 



(k 2 



. 
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Using formulations above we are able to calculate 
the primordial power spectra with mode by mode 
integrations [l4l 134 l35| . Regarding the choices of the 
model parameters: initial values of <pi and 02, mi and 
ni2, we notice that 0i is arbitrary with a weak prior to 
provide enough number of e— folding to solve the flatness 
problem|24l l2fi|. 02 mainly determines which Pr corre- 
spond to the cosmological scale and the ratio of mi to 
77i2 determines the shape of Pr with the absolute value 
of 777,2 fixed by the WMAP normalization. For different 
values of 02 , the corresponding number of e — folding to 
CMB scales will differ and the shape of Pr will also get 
changed, as in the case of the single field inflation. The 
amplitude of Pr(/s) is to be determined by observations. 
In Fig^ we show initial power spectra as a function of 
\n(k/kf). In the numerical calculation |3(j we have set 
mi = 87712 and 02 = 3.3A/p; at the onset of inflation. 
This gives rise to N(kf) — 59.6. We find such a value of 
N(kf) is acceptable for fitting to WMAP data below. We 
also show in Fig^ the behavior of the slow rolling(SR) 
parameters defined by e = —H/H 2 and 5 = a/Ha . 
One can see that these parameters change dramatically 
- this is why we use numerical calculations instead of the 
Stewart-Lyth analytical formula |33L l35l l3^| . 

Now we fit the WMAP data with the primordial spec- 
tra in Fig^ Our modified version of the publicly avail- 
able CMBFAST H3 is based on Version 4.2 [U and we 
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FIG. 2: CMB anisotropies with primordial spectra shown in 
Fig0- The error bars are taken the same as Ref.0- The up- 
per and lower dashed lines show the 1-a confidence levels for 
lognormal distributions with cosmic variance limits. The mid- 
dle solid lines show the models with the lowest quadrupoles. 



FIG. 1: Initial power spectra and slow rolling parameters as 
a function of ln(fe/fe/). N{k f ) = 59.6 



have used the "HP" choice to give exact CMB TT and 
TE power spectra. We also run with CAMBpoL HP] for 
a crosscheck on our results. We use a similar method to 
Ref . pill . and set J7a and lnfc/ as free parameters in our 
fit. Denoting k c = 7.0 x 70./3/10 5 w 1.6 x 10" 3 Mpc~\ 
we use 101 and 251 grid points with ranges [0.68,0.77], 
and [—19.6,5.4] respectively for Q\ and ln(fc//fc c ). At 
each point in the grid we use subroutines derived from 
those made available by the WMAP team to evaluate 
the log likelihood with respect to the WMAP TT and 
TE data |3|. Other parameters are fixed at Qhh 2 — 0.022, 
il m h 2 = 0.135 and r c = 0.17 and il to t = 1 @ The over- 
all amplitude of the primordial perturbations has been 
used as a continuous parameter. Differing from Ref.|24T] 
and Ref.[2{|, we have included the tensor contributions 
in our fit to CMB. And for comparison, we also run the 
code without tensor. 

In Fig. [5] we show the resulting best-fit models ob- 
tained from the grids for the model considered. In the 
plots we take the same error bars on the binned WMAP 
results as Ref.0. The regions between the two dashed 
lines are given by 1-a confidence levels for lognormal dis- 
tributions as well as the cosmic variance limits. The cos- 
mic variance factor is 1 ± ^2/(21 + 1). The middle solid 
lines show the models with the lowest quadrupoles. We 
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FIG. 3: Two-dimesional contours in the ln(fc//fc c )- CIa plane 
for our grids of model. k c ~ 1.6 x 10 -3 Mpc -1 . The regions 
of different color show 68.3%, 95% and 99.7% confidence re- 
spectively. 



get the minimum \ 2 = 1429.2 when not including tensor 
and a slightly larger x 2 — 1429.7 with tensor for the best 
fit values. One can see the CMB low multipoles do have 
been suppressed in our model. 

In Fig. |3| we plot the resulting x 2 values as functions of 
and ln(A://fc c ). The contours shown are for A% 2 val- 
ues giving one, about two, and three a contours for two 
parameter Gaussian distributions. We find that the pri- 
mordial spectrum with a feature is favored at more than 
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FIG. 4: One-dimesional marginalized distributions for the val- 
ues of ln(fc//fc c ). 

3-cr level. However, when neglecting the tensor contri- 
bution, the significance reduces to only about 2-a level. 
In our model the ratio of tensor to scalar r = 8P g /Pn 
reaches its minimum at N ~ N(kf) and it gradually in- 
creases for smaller TV (larger k) and grows rapidly for 
larger N (smaller k). And our N(kf) corresponds to the 
region with the lowest value of r, which is consistent with 
WMAP group's analysis that large tensor contribution is 
disfavored by current CMB observations Q|. The 1-a 
regions in Fig. 2] and Fig. [21 differ slightly when with and 
without tensor since r is around its minimum in both 
cases. 

We marginalize over J7a to obtain the one-dimensional 
probability distributions in ln(kf/k c ) shown in Fig. 0] 
For the spectra in Fig^ when neglecting tensor contri- 
butions, we get hi(kf/k c ) = 0.8 with the maximum like- 
lihood, corresponding to kf = 3.6 x 10~ 3 Mpc -1 . We 
also have kf ~ at 2a level. When taking into account 
tensor contributions, the maximum likelihood value of 
ln(kf/k c ) shift to 1.3 and we get ln(fc//fc c ) = 1.3^q' 4 
at lcr, 1.3j^»j at 2a and 1.3j1;'o a ^ 3c, corresponding 
to k f = 6.01 3 ,;? x lO- 3 Mpc- 1 , 6.0±fj x lO^Mpc" 1 
and 6.OI5 2 x 10 _3 Mpc _1 respectively. The difference 
in x 2 between the peak in the distributions and at 



ln(kf/k c ) = —19.6 is found to be A% 2 = 6.1 when not 
including tensor and A\ 2 — 18.3 when with tensor. 

When considering tensor contributions in the two di- 
mensional contour between ln(kf/k c ) and the normal- 
ized factor of the primordial scalar spectrum at k = 
0.05 Mpc^ 1 , we obtain h\(kf/k c ) = -0.8 ~ 1.9 and 
Pr(0.05/M P c) = 2.48 - 2.53 x 10~ 9 at 2a level and 
we get mi = 1.3 ~ 1.4 x 10 14 GeV. In contrast to the 
power law primordial spectrum with constant ns we run 
a similar code: we fix fli,h 2 = 0.022, il m h 2 = 0.135 and 
t c = 0.17 and f2 tot = 1, varying tt\ and ns with ranges 
[0.68,0.77], [0.91,1.07] and get a minimum X 2 = 1432.7. 
To characterize the pure power law primordial spec- 
trum, one considers two parameters: ns and the am- 
plitude. For our double inflation model four parame- 
ters are introduced to give the exact scalar and tensor 
spectra: m 1; mi/m2, N(kf) (or equivalently (f>2 at the 
onset of inflation) and ln(fcy/fc c ). This indicates our 
double inflation model is favored at ~ 1.2cr compared 
with power law ACDM model^l- I n general, primor- 
dial power spectra with a feature or cutoff do generate 
a lower CMB TT quadrupole (which, however may not 
be sufficient) 0, 0> |25[ . A cutoff primordial spectrum 
also as pointed out in Ref . [2{| makes the CMB TE multi- 
poles lower. When combining these two effects, however, 
our calculations show that the primordial spectrum with 
a feature can work better than models with power law 
primordial spectra. 

In conclusion, we have studied the possibility of sup- 
pressing the low multipoles in the CMB anisotropy with 
a model of double inflation. Our results show that with 
mi ~ 10 14 GeV and m 2 - 10 13 GeV which lies in the 
parameter space required by neutrino physics in the sce- 
nario of sneutrino innation|2jj|30(, this model fits to the 
WMAP data better than the standard power-law ACDM 
model. 
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